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*> 1 Introduction 

§ . The CalcPHEP collaboration joins the efforts of several groups of theorists known very well 
CN ! in the field of theoretical support of various experiments in HEP, particularly at SLAC and 



LEP, (see, for instance [|I|], |2| and ||). The first phase of the CalcPHEP system was realized 



in the site \http://brg.jinr. ru/\ in 2000-2001. It is written mostly in F0RM3, |4j]. In this talk, 



we will describe the present status and our plans for the realization of next phases of the 
CalcPHEP project aimed at the theoretical support of experiments at modern and future 
Oh! accelerators: TEVATRON, LHC, electron Linear Colliders (LC's) i.e. TESLA, NLC, 
CLIC, and muon factories. Within this project, we are creating a four-level computer 
system which eventually must automatically calculate pseudo- and realistic observables 
for more and more complicated processes of elementary particle interactions, using the 
principle of knowledge storing. Upon completion of the second phase of the project, 
started January 2002 with duration of about three years, we plan to have a complete 
set of computer codes, accessible via an Internet-based environment and realizing the 
complete chain of calculations "from the Lagrangian to the realistic distributions" at the 
one-loop level precision including all 1 — > 2 decays, 2 — *> 2 processes and certain classes of 
2^3 processes. 



1.1 CalcPHEP group 

The CalcPHEP group was formed in 2001 in sector N°l NEOVP LJAP. 

During the first phase of the project in 2000-2001, the CalcPHEP group created the 

site brg . j inr . ru, where the development in two strategic directions is foreseen: 

1. Creation of a softwear product, capable to compute HEP observables with one-loop 

precision for complicated processes of elementary particle interactions, using the 
principle of knowledge storing. Application: LHC. 

2. Works towards two- loop precision level control of simple processes: 1 — > 2, 1 — > 3 
and 2^2. Application: GigaZ option of electron LC's. 
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1.2 A little bit of history 



There are two historical sources of CalcPHEP project: 

1. From one side it roots back to many codes written by Dubna group aimed at a theo- 
retical support of HEP experiments in the past: 

1975 - 1986: support of CERN DIS experiments (BCDMS, EMC, NMC), creation 
of program TERAD; support of CERN neutrino experiments (CHARM-I, CDHSW and 
CHARM-II), creation of programs NUDIS, INVMUD, NUFITTER. 

1983 — 1989: Foundation of the DZRCG - "Dubna-Zeuthen Radiative Correction 
Group", creation of EW library DIZET; creation of the program ZBIZON — the fore-runner 
of ZFITTER 0. 

1989 — 1997: support of the DIS experiments at HERA, creation of the program HECTOR; 
participation in SMC experiment at CERN with the program fiela. 
1989 - 2001: Theoretical support of experiments at LEP, SLC (DELPHI, L3, ALEPH, 
OPAL and SLD). 

2. From the other side, a monograph "The Standard Model in the Making" was writ- 
ten |J. While working on the book, the authors wrote hundreds of "book-supporting" 
form-codes, which comprised the proto-type of future CalcPHEP system. 

Like well known codes of LEP era: T0PAZ0 [0, ZFITTER [§, KXMC [§, CalcPHEP is 
supposed to be a tool for precision calculations of pseudo- and realistic observables. Let's 
remind these definitions that arose in depth of LEP community: 

Definition 1 Realistic Observables are the (differential) cross-sections (more general event 
distributions) for a reaction, e.g. 

e+e~ -> (7, Z) -> fj(n~/) 

calculated with all available in the literature higher order corrections (QCD, EW), in- 
cluding real and virtual QED photonic corrections, possibly accounting for kinematical 
cuts. 

Definition 2 Pseudo-Observables are related to measured quantities by some de-convo- 
lution or unfolding procedure (e.g. undressing of QED corrections). The concept itself of 
pseudo-observability is rather difficult to define. One way say that the experiments measure 
some primordial distributions which are then reduced to secondary quantities under some 
set of specific assumptions (definitions). 

Z decay partial width represents typical example of pseudo-observables, i.e. it has to 
be defined. At the tree level, we define it as a quantity described by the square of one 
diagram: 
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2 LEP, Precision High Energy Physics and its Future 



One may say that during recent years a new physical discipline was born. We call it 
PHEP, Precision High Energy Physics. Experimentally, it finally shaped in the result of 
glorious 12 year LEP era: measurements at Z resonance in 1989 — 1995, and reaching 
an unprecedented experimental accuracy < 10~ 3 , and measurements above Z resonance 
in 1995 — 2000, at higher energies, where high enough experimental accuracy was also 
reached < 1%. By 2/11/2000 LEP2 possibly saw hints of "God blessed" particle — Higgs 
boson, but was stopped, unfortunately, mainly due to lack of financing. 

For the first time huge HEP facility challenged for theoreticians to perform calculations 
with uncertainty better than experimental errors of O (10~ 3 ) and, eventually, efforts of 
many groups of theoreticians allowed the achievement of the theoretical precision of the 
order 2.5 • 10~ 4 at the Z resonance and 2 — 3-10 4 at LEP2 energies. 

This, in turn, greatly contributed to the success of precision tests of the SM, the main 
result of LEP era, which laid the foundation of the Precision High Energy Physics. This 
is why our project got this suffix PHEP. 

2.1 Future of PHEP 

PHEP has good perspectives and after the end of LEP. Several Input parameters of the 
Standard Model (SM) are expected to be improved in near future. 
Recent discrepancy in the muon amm: 

116591661(114) x 10~ n 
116592023(151) 
362(189) 2a difference 

exp. error (151) should be improved soon up to ~ (50), (1) 

necessitates an improvement of the knowledge of the hadronic contribution of Aa^ (M 2 ^ 
to the running e.m. coupling. An experimental input for a (e + e~ — > hadrons) at cms 
energies (1-4 GeV) is expected from BES-II, BEPC (Beijing), VEPP2000 (Novossibirsk) 
and DAFNE at cms energies around 0-meson. 

Very important should be projected improvements of mass measurements: M w , m t . 
LEP1 finished with indirect result for the top mass: m t = 169tg° GeV; while LEP1 © 
TEVATRON constraint yielded: m t = 174.5^2 GeV. 

LEP2 reached for W mass M w : M w = 80.450 ± 0.039, in the direct measurements 
and M w = 80.373 ± 0.023 as indirect result. 

TEVATRON in RUN-I reached: M w = 80.454 ± 0.060 GeV, m t = 174.3 ± 5.1 GeV. 
Much better precision tags are expected to be reached at TEVATRON, RUN-II (re- 
cently started): AM W ~ 20MeV, Am t ~ 2 GeV; and later at LHC (not so sooner than 
in 2006, however): AM W ~ 15 MeV, Am t ~ 1 GeV. 

Where, when and with which mass Higgs boson might be discovered? 

— TEVATRON has a serious chance to see Higgs up to mass 180 GeV; however it will 
require very high integrated luminosity: f £ > 5/6 -1 ; 

— LHC, will cover all allowed mass range up to 500 GeV (not so soon, after 2007); 

— LC's and muon factories (after 2010-2012). 



a E il xp (Average) 
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New horizon of PHEP will be opened with experiments at electron LC's: TESLA 
(DESY) particularly with GigaZ option, i.e. coming back to Z resonance with statistics 
10 9 ; CLIC (CERN); JLC (KEK), NLC(SLAC, LNBL, LLNL, FNAL) and Muon Storage 
Rings (Higgs Factory) — all that more than in ten years from now. 

One expects fantastic precision tags there in: 

- Asin 2 # e// ~ 0.00002; 

- AM W ~ 6 MeV, Am t ~ 100 - 200MeV; 

- AM H ~ 100 MeV (from e+e~ -> ZH); 

- and detail study of Higgs boson properties. 

Given our LEP1 experience one should definitely state that 2-loop precision level 
control will be absolutely necessary for the analysis of these data! 

One may conclude that PHEP has a bright future: all future colliders — TEVATRON, 
LHC, electron LC's (TESLA, NLC, CLIC) and muon factories will be, actually, PHEP 
facilities! For data analysis, they will surely require qualitatively new level of both theo- 
retical predictions and principally new computer codes. 

3 Necessary notion 

In order to understand the language of CalcPHEP one has to introduce many notions and 
notations. 

3.1 Input Parameter Set, IPS 

The Minimal Standard Model (MSM), contains large number of Input Parameters: 
25 = 2 interaction constants a and a s 

© 8 mixing angles (CKM and possible lepton analogs) 

© 15 masses (12 fundamental fermions and 3 fundamental bosons Z, W, H). 
However, the number 25 is minimal. MSM is unable to compute its IPS from first 
principles; MSM is able to compute any observable 0^ xp in terms of its IPS: 

Of p (measured) <-> 0* heor (calculated, as a function of IPS) . (2) 

This is the way how precision measurements set constraints on IPS. 

3.1.1 Number of free parameters in fits of Z resonance observables 

At Z resonance, not all 25 parameters matter. Actually only 5 parameters: 

A«i 5) (Mj), a s (M 2 z ), m t , M z , M H , (3) 

which we call the Standard LEP1 IPS, matter. 

Using M z , measured at Z peak itself with the precision ~ 2 x 10 -5 , and also reach 
information from the other measurements for: 

a s (M 2 z ), m t , M w , (4) 

we approach one-parameter fit, with Higgs boson mass M H being the only fitted parameter. 
The result of such a fit was shown in the Blue band figure, the most celebrated LEP era 
figure, derived with the aid of T0PAZ0 [0 and ZFITTER [|] codes. 
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3.2 Quantum Fields of the SM 

Here we sketch all fundamental quantum fields of the SM in one of the most general 
gauges — with three arbitrary gauge parameters £ A , £ z , £. 
Three generation of fermions or matter fields: 



/ = 
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D 
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possess masses, m/, charges, Qf, and third projections of weak isospin, If 
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Gauge fields: 



Vector bosons 
WWW A 

/ W^V Z{M Z ) 
/ WW W ± (M W ) 
Gluon 

possesses strong interaction 
WW g 



Unphysical scalars Faddeev-Popov ghosts 

► Y' 



- 0° 

- ± 



yz 

X ± 



possess physical charges and physical masses 

possess physical charges and unphysical masses 

and unphysical charges. 

Higgs field: 

H (M H ) is a scalar, neutral, massive field. 



3.2.1 The Lagrangian in gauge, Feynman Rules 

At the ground level of CalcPHEP system one has this Lagrangian 

C = £(IPS of 25 parameters, 17 fields, 3 gauge parameters), 
from which one derives primary Feynman rules for vertices. 

5 



(5) 



3.2.2 Propagators in gauge 

Here we list propagators in gauge, the other important bricks of CalcPHEP system. 
Propagator of a fermion, /: 



— ij> + rrif 

f p 2 + rr?t 



Vector boson propagators: 



V 



A WWW ^{^ + (C-^) P -fj 



p 2 + M 2 ^ z I p 2 + i\Ml 

p 2 + M 2 \ 5flu ' V* V p 2 + i 2 M 2 

Propagators of unphysical fields: ^ C 

Y A p 

1 . £ 



00 p 2 + £2 M 2> yz p 2 +£|M 2 

1 f 

"V"" p 2 +f 2 ^ ' ""' i±"""'"" p 2 + e 2 ^ 

Propagator of the physical scalar field, if-boson 

1 



H p 2 + M 2 H 

3.3 Scalar Aq, Bq, etc functions 

For calculation of one- loop integrals CalcPHEP uses the standard scalar A , B , Cq and 
D functions |J. 

One-point integrals or A functions, are met in tadpoles diagrams: 




We give its defining expression: 



zvr 2 A (m) = / d n q \ . . (G) 



gZ _|_ m Z _ ^ 

and the answer in the dimensional regularization 



A (m) = m 2 {--_ - 1 + ln^P) + O (e) , (7) 



r /U 2 



where the ultraviolet pole is: 

1 2 

- = 7 — ln7r, n = 4 — e. (8) 

Two-point integrals or Bo-functions are met in self-energy diagrams: 




m 2 

We limit ourselves by giving its defining expression: 

{ 1 

m 2 B (pl;m 1 ,m 2 ) = J rfn g^- > 

d = q 2 + ml - ie, d\ = (q + pi) 2 + m\ - ie . (9) 
Three-point integrals, C functions, are met in vertices: 




Its defining expression reads: 

{ 1 

m 2 Co (p 2 1 ,p 2 2 ,Q 2 ;m l} m 2} m 3 ) = J rfn g^-^ , (10) 

do = q 2 + m\ - ie, d x = (q + pi) 2 + m\ - ie, d 2 = (q + p x + p 2 ) 2 + m\ - ie, (11) 

where Q 2 = (p\ + p 2 ) 2 is one of the Mandelstamm variables: s, t or u. 
Four-point integrals, D -functions, are met in boxes. 

Presently, CalcPHEP knows ALL about reduction of up to four-point functions up to 
third rank tensors and of the so-called special functions, which are due to peculiar form 
of the photonic propagator in the gauge, see 0. 

3.4 Processes in the SM 

One should be aware of a hierarchical classification of processes accepted in CalcPHEP 
and of a relevant notion of independent structures, or independent amplitude form factors, 
which number is deeply related to the number of independent helicity amplitudes by which 
a process may be described (below we present these numbers for unpolarized cases). 
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3.4.1 Decays 1 -> 2 

There are B — > // and 3-B decays: 



• 


if - 




(one structure) 


• 


z - 


- 77, (7 - //) 


(three structures) 


• 


w - 


- //', (f - w + b) 


(four structures) 


• 


H - 


-> zz, w + w~ 




• 


Z - 


* w+w~ 





3.4.2 Processes 2 2 

There are 2f — > 2f processes, which in turn are subdivided into Neutral Current (NC) 
and Charged Current (CC) ones: 

• NC: // — >• (7, Z, H) — > /'/' (4,6) 10 structures depending on whether initial and 
final state fermion masses are ignored 

. CC: fj 2 - (W) -> fj 4 

Next, there are many processes of a kind Vf — > f'V, in particular 

• compton-effect: — > je, Z — > //7 

• e+e- -> W^+W^-, ZZ, Z 7 ,77 

Decays 1 — >• 3 are cross-channels of the previous processes and their one-loop descrip- 
tion in terms of independent objects, mentioned above, one gets for free. Present level of 
CalcPHEP has a lot of preparations for all above processes, but far not all is put into the 
working areas of the site brg . j inr . ru. 



3.4.3 Processes 2^3 

They comprise a very reach family, for instance: 

• e+e--(7,Z,ff)-/77- 

Their implementation is one of main goals of the second phase of CalcPHEP project. 
Corresponding decays 1 — > 4 are again cross-channels of the previous processes and need 
not be studied separately. 



3.4.4 Processes 2^4 

To this family belongs 4 fermion processes of LEP2. Their study is not foreseen at the 
second phase of CalcPHEP project, but might be a subject of its third phase. 
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4 Building Blocks and knowledge storing 

4.1 Simplest decay: Z —* ff 

4.1.1 Amplitude of Z — > // decay at tree level 

Its tree level diagram was already presented at the end of Section 1.2; the corresponding 
amplitude reads: 



v zf7 = (2 



,4 
7T Z 



w 

2c„, 



■7/, 



/f(i + 75 )-2g /S 



with vector and axial coupling constants: Vf 

of the two structures in Eq. (|12"D, which might be termed as L and Q structures, corre- 



, a f 



'12) 



1^. Note appearance 



spondingly. Note also, that Eq. (12) as well as all below, are written in Pauli metrics that 
is used by CalcPHEP. 



4.1.2 Amplitude of Z — ► // decay with loop corrections 

It might be schematically depicted as a sum of one-loop vertices and counter terms: 

7/pi 7/ 7, 

z ^ z ^ z 

+ 

f\P2 f\ f 

In the most general case (but for unpolarized study) the one-loop amplitude may be 
parameterized by the three scalar form factors: 

,3 






9" 



Vf// = (2vr)S ; 

V ; 16tt 2 2c„ 



■7m 



ilf ] F Ll+ - 2tQ f s 2 w F Q + m f ( Pl - p 2 ) F D 



(13) 



Given similarity of Eqs. (0) and (|13|) , the latter is called sometimes Improved Born 
Approximation (IBA) amplitude. 

4.1.3 QED diagrams and corrections 

The QED diagrams comprise gauge invariant subsets, this is why they are considered 
sometimes separately: 

7/ 7/ 1, 




+ 




+ 




A A / 

Their contribution to the partial Z widths, in the case when no photon cuts are imposed, 
reads: 



-pQED 
1 / 



(14) 
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4.2 Process e + e 



Coming to a more complicated case of a 2f — > 2/ process, we will illustrate how building 
blocks, derived for a study of a lower level process, might be use at a higher level. 

4.2.1 Tree-level diagrams and amplitudes of e + e~ — > // 

Consider first the two tree-level diagrams with 7 and Z exchanges in order to introduce 
basis of relevant structures. 



+ 





A 



Born 



e 2 Q e Qf 



A 



Born 



7m ®7^ 



4s 2 c 2 

14/ IV 



7 A1 (u e + ^75) ® 7m( u / + a /7s) 



4s 2 c 

W W 



7„ 



(15) 



where 7± = 1 ± 75 and symbol <S> stands for a short-hand notation 

7^ («i + «i75) ® 7^ («2 + 0275) = v (p+) 7 M (ui + 0175) u (p_) it 7„ (t> 2 + a 2 7s) u (<7+) 
and 



Xz{s) 



(16) 



s-M 2 +isr z /M z ' 
This amplitude is characterized by four structures: 

LL = 7^7+ ® 7^7+ ) ^Q = 7 m 7+®7m> <2 L = 7 M ® 7^7+ , <2Q = 7 M ® 7 M • ( 17 ) 

4.2.2 One-loop amplitude for e + e — > // 

"Dressed" with one-loop vertices and counterterms, the 7 and Z exchanges may be sym- 
bolically depicted as: 



+ 
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And similarly for the initial state vertex: 






+ 



Where one can easily recognize building blocks already known from the calculation of 
one-loop radiative corrections for Z decay, however, now we need to dress 7 — > ff vertex 
too and add into consideration "dressing" of propagators: 



+ 



To complete calculations of one-loop EWRC for the process e + e 
add WW and ZZ boxes: 

e+ W d 
« *AAAA/\AAA/t « 






f f one should 



u 



+ 



— ► — *A/VW\AA/V»— »— 

e" w d 
e+ (Z,j) / 
--•-WWWWWIM- 




e u 
e + (Z, 7 ) 7 



/ 



+ 



<W\A/\AA/Ww — *— 

(z,j) f 




e- (Z, 7) / 

ZZ, 77 and Z7 boxes comprise gauge invariant subset of diagrams. Moreover, 77 and 
Z'-f boxes QED, vertices and QED bremsstrahlung for NC 2f — > 2f processes often are 
separated into a gauge invariant QED subset of diagrams. 

Virtual QED one-loop diagrams together with four QED bremsstrahlung diagrams 
form an Infra-Red Divergence (IRD) free subset. 

This example clearly shows how the principle of knowledge storing is implemented 
within CalcPHEP project: one starts from the simplest decays and collects all relevant 
building blocks, BB's (off-shell with respect to boson mass). Then one moves to next 
level of complexity where all BB's computed at the previous level are requested, but on 
top one needs more complicated objects (here boxes). 

This strategy was realized in our recent calculations of the EWRC to the e + e~ — * ff 
process, which are completely done with the aid of CalcPHEP system 0. There is another 
study accomplished with CalcPHEP [ 10 1 . 
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5 Status of the project 



Before discussing what is already available at the site brg.jinr.ru, we present some 
general information about CalcPHEP system. 



5.1 Basic information about CalcPHEP, keywords 

CalcPHEP is four-level computer system for automatic calculation of pseudo- and 
realistic observables (decay rates, event distributions) for more and more complicated 
processes of elementary particle interactions, using the principle of knowledge storing. 
At each of the four levels there are: 

1. Codes (written in FORM3), realizing full chain of analytic calculations from the SM 
Lagrangian C SM to the Ultra Violet Free Amplitudes, UVFA, parameterized by a 
minimal set of scalar form factors; 

2. Codes (written in FORM3), realizing analytic calculations of a minimal subset of 
Helicity Amplitudes, HA's, followed by an automatic procedure of generation of 
codes for numerical calculations of HA's (presently FORTRAN codes, and in a near 
perspective C ++ codes). 

3. Codes, realizing the so-called "infrared rearrangement" of HA's. This is needed if the 
multiple photon emission is being exponentiated at the amplitude level. Currently, 
bremsstrahlung photons are added in the lowest order and the third level is skipped. 

4. Codes, that use HA's derived at the second (or third) level together with tree-level 
HA's for one-photon (or multiple-photon) emission, within a Monte Carlo event 
generator, which is supposed to compute realistic distributions (presently FORTRAN 
codes, and in a near perspective C ++ codes.) 

It is an Internet based and Database based system. The latter means that there 
is a storage of source codes written in different languages, which talk to each other. They 
are placed into a homogeneous environment written in JAVA. 

It follows Intermediate access principle i.e. full chain "from the Lagrangian to 
realistic distribution" should work out completely in real time, if someone requests this, 
however, it is supposed to have several "entries", say after each level, or just providing 
the user with its final product — a Monte Carlo event generator. 



5.2 Some technical data about CalcPHEP 



1 . Address \http://brg. jinr. ru/\ 



2. For realization of the site one used: 

— Apache web server under Linux, 

— form3 compiler, 

— mySQL server for relational databases. 

3. In the current version, user-interface is realized with the use of PHP. 

4. Nowadays, everything is being rewritten in JAVA in order to reach better "interac- 
tivity" and to use reach possibilities of already written in this language libraries. 
Main goal of this rewriting is to create a homogeneous environment both 
for accessing our codes from the database and for offering a possibility for 
simultaneous work of several members of the group and external users. 
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5.3 Present and nearest versions of CalcPHEP system 

In 2001, we released two test-versions of CalcPHEP: 

1. vO.Ol from March'01 realizes analytic calculations of one-loop UVFA for de- 
cays 1 — > 2 ( level- 1). 

2. v0.02d from March'02 returns numbers for one-loop decay widths (levels-1,2) via 
temporary bypass of level 4. It realizes also levels-1,2 for 2f — > 2f NC process. 

3. One has very many almost finished "preparations" for the other processes 2^2 
and decays 1 — ► 3 (level-1). All this should comprise v0.03 of Summer 2002. 

4. An active work is being realized on implementation of level-4 for decays 1 — > 2, this 
should complete full chain "from the SM Lagrangian to pseudo-observables" for the 
simplest decays. 

5. There are many problems to be solved at the second or later phases of the project. 
Among them one should mention: 

- automatic generation of Feynman Rules from a Lagrangian, 

- automatic generation of topologies of Feynman diagrams, 

- graphical representation of the results. 



6 Conclusion 

At a Symposium in honor of Professor Alberto Sirlin's 70th Birthday was said: A new 
frontier is as the horizon: most likely it is goodbye to the one man show. Running a new 
Radiative Correction project will be a little like running an experiment JI1|/ . 

Indeed, projects of such a kind as CalcPHEP are definitely long term projects. Remem- 
ber, that ZFITTER took about 12 years, about the same time exists already FeynArts [12[ . 



Our nearest goal is the realization of the second phase of the project upon completion 
of which we plan to have a complete software product, accessible via an Internet-based 
environment, and realizing the chain of calculations "from the Lagrangian to the realistic 
distributions" at the one-loop level precision including some processes 2 — > 3 and decays 
1 — > 4. Plans also assume to perform an R&,D for the third phase of the project (see 



also |L3] and successive developments of [|T4|-|0) which should begin in 2004. 



Second phase is basically oriented on a common work of theoreticians of the Dubna 
group and the Knoxville-Krakow collaboration [0] . 

United group proposes to realize in 2002-2004 an important phase of CalcPHEP 
project: oriented toward a merger of analytic results to be produced by Dubna team 
with MC event generators to be developed by Knoxville-Krakow collaboration^. 

Among most important milestones of first year, one should mention: realization of the 
levels 2-4 for the simplest Z(H, W) — > // decays; completion of level 1 for the radiative 
Z decay, Z — > ffj, work on which is already under way; completion of levels 2-4 for the 
radiative Z decay. 



x In this connection it is necessary to emphasize that any future code aimed at a com- 
parison of experimental data with theory predictions should be a MC generator, since the 
processes at very high energies will have multi-particle final states that make impossible 
a semi-analytic approach used at LEP within ZFITTER project. 
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